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Abstract 



The recent proposals concerning the usage of the real polarized gamma beam, ob- 
tained by the Compton backscattering of the laser photons off the electron beams 
from either the linear or circular accelerators were considered. The heavy quark 
photoproduction process giving a unique opportunity to measure polarized gluon 

00 ■ distribution was investigated. 
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^ ■ 1 Introduction 



A spin crisis arose when the first determination of 

AS = Aw + Ad + As 
was found to be much smaller than expected [1], where 



Af = J dxAf(x,Q 2 



and f(x,Q 2 ) are the polarized quark spin distribution functions. The recent 
value of the world average for AE is approximately 0.3 ± 0.06. It is much 
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smaller than the relativistic quark model prediction of 0.6 [2]. Among the 
number of explanations of the EMC results on the longitudinally polarized 
proton structure function, that of the perturbative QCD approach occupies a 
prominent position and opens a new domain of tests [3]. It has been shown 
that in perturbative QCD predictions of the quark-parton model concerning 
the singlet axial charge contribution to If need to be modified because of 
the 75 anomaly of the flavour singlet current J^. In this approach the main 
question is the size of the polarized gluon distribution. But there is a great 
uncertainty in theoretical estimates on the magnitude and re-dependence of 
AG(x), especially in the moderate x-region x > 0.1. One calculation, per- 
formed in the framework of MIT Bag Model, predicts negative AG(x) value, 
thus even sharpening the problem [4]. 

In order to measure gluon contribution to the nucleon spin, we must select 
a process where the familiar lowest order graphs of deep inelastic scattering 
from a single quark are suppressed. Analysis of the experimental data shows 
that there is no significant charm content in the nucleons [5]. Then one may 
hope to determine polarized gluon effects in the process of the charmed quark 
photoproduction. Due to the comparatively large mass of the c-quark the 
leading mechanism for obtaining a cc- pair in the final state is the hard photon- 
gluon fusion (PGF) process 



7 + #->c + c (1) 

and produced cc- pair may form J/ip meson or fragment separately into open 
charm states. 

A rich spin physics programme was proposed by COMPASS collaboration at 
CERN [6]. A large part of it is devoted to measurement of gluon polariza- 
tion by tagging D°, D* + -mesons produced in the 7*iV collisions. Produced 
D°, _D* + -mesons will be reconstructed from their two-and three-body de- 
cays into hadrons. It permits to impose constraints on the invariant masses 
of K 7i 7r, K 7r, .. subsystems, which effectively rejects a background, espe- 
cially in the case of D* + -tagging. A statistical precision on the polarization 
asymmetry measurement is expected to be about 0.05. 

A proposal has been made [7] to determine AG(x) from the measurement 
of the asymmetry in the charm photoproduction process in the scattering of 
polarized real photons off polarized fixed target. In this paper we consider a 
possibility to realize this experiment at the TESLA machine [8], and discuss 
briefly an opportunity to perform it at SLAC and LEP2. 
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2 The Real Gamma Beam 



2.1 Linear Accelerator 



The scheme of the proposed experiment looks as follows. Circularly polar- 
ized laser beam with photon energy c^o = 3.3 eV ( Cul5 laser) is scattered 
off the 250 (50) GeK-electrons provided by TESLA (SLAC) [9]. Throughout 
this subsection the numbers in parentheses refer to the SLAC gamma beam 
parameters. The Compton back-scattered photon beam has a spectrum 7(2/) 



/ , 1 dN Ar r dN 
1{V) = N X ^ N = J dy ^ 
dN 1 . 

— = - + l- y -4r 1-r , 2 

dy l-y 

and follows closely to the trajectory of the primary electron beam. Here the 
ratio of the hard photon energy uj to the electron energy E e ranges as 



y = ^r, 0<y<y max = -^—-0.927 (0.717), 

k = ~ 12.638 (2.528), r = — - — -. (3) 

mj K{l-y) 

Had the laser beam been circularly polarized, the 7-beam would be polarized 
too and the 7-beam helicity is given by 



^ JV 7 (y)' 
B(y) = -X ph (2r-l)(^ + l-y). 

y 

where X p h denotes the laser photon helicity. Due to the unique dependence of 
both 7-beam energy and polarization on the scattering angle # 7 between the 
incident and scattered photon 



6I»~— -(« + !), 

it is possible to obtain almost monochromatic gamma beam with energy 
0.99u max < uj < uj m ax and polarization nearly equal to unit by selecting pho- 
tons with # 7 < 0.759 (1.93) • 10~ e rad. Taking the distance between the con- 
version region and collimator 100m, we obtain that this angle corresponds to 
the diameter of the selecting slit d = 152 (386)//m. 
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In our further estimation on the luminosity of gamma beam scattering off 
a polarized fixed target, we shall follow to the method outlined in [10]. We 
determine the optimal number of converted photons iV 7 by the requirement 
to obtain one event per collision 



/3iV 7 T n a 7P = 1, (4) 

where f3 is the fraction of the photons coming through the slit, T n stands for 
the target nucleon density and <t 7P is the total cross section of gamma-proton 
collisions. In the case of 1% 7-beam monochromaticity (3 ~ 4.22 (2.02) • 10 -2 , 
for deuterated butanol target with the length about 40 cm the density is T n = 
4-10 25 cm- 2 and a rp ~ lOOufr at high energies, so we get iV 7 = 5.92 (12.36) -10 3 . 
Because the number of electrons in a bunch N e = 5.15 (3.5) • 10 10 for TESLA 
(SLAC), we obtain the necessary conversion coefficient 



# = ^ = 1.15 (3.53) -10- 7 . 

1 ' e 



The number of photons required to provide the scattering of each electron 
with the laser photon is defined from 

^ = 1, (5) 

where no is the number of photons in a laser pulse, a c — 1.1 (2.59) ■ 10~ 25 cm 2 
is the total Compton cross section and S e ff is the effective area of the photon 
and electron beams intercept. Choosing S e ff = 4 • 10™ 6 cm 2 , one can easily get 
the laser energy per pulse 

A = n uj ~ 19.2 (8.152) J, 

corresponding to the total electron conversion. Since the optimal number of 
converted photons iV 7 , defined by (4), is much less than N e , we need not the 
total conversion of electrons and the required pulse energy is suppressed by 
the factor of K 

A = K A ~ 2.21 (2.88) • 10~ 6 J, 
which is accessible with modern laser technology. 

Finally we shall discuss a choice of the laser pulse frequency. It should evi- 
dently coincide with the linac frequency f pu i se when number of the bunches 
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in the electron beam is equal to unity. If the electron beam has a multybunch 
structure, one may either choose the repetition rate of laser pulses to be equal 
to f p U i S e • n b or use a mirror system [11], which evidently decreases the laser 
pulse frequency. Since the detailed discussion of an experimental setup is be- 
yond the scope of our paper, we choose for simplicity the laser pulse frequency 
equal to f pu i se ■ n b . 

For the linear accelerators the integrated luminosity per a year of operation 
takes the form 

KL = 10 7 • frep f3N e K T„, (6) 
frep fpuise ■ n b , 

where f rep is the collision frequency. Taking into account the Eq.(4), one can 
estimate L l ^ ac as follows 



= 10 7 -^~ 0.8 (0.012) fb-\ (7) 

(J^p 

2.2 Circular Accelerator 

Now we briefly consider the possibility to obtain the real gamma beam at a 
circular electron accelerator, say at LEP2. All calculations are done in analogy 
with the TESLA case, except a few details. First of all, f rep now takes the form 



frep = fpuise • U b = — — ■ U b = 44 980 Hz, 

lixti 

where the number of electron bunches n b — 4, c is the speed of light, R is 
the ring radius and f pu ise is the circulation frequency in this case. It is evident 
that for the circular accelerator f rep should coincide with frequency of the laser 
pulses. The next remark concerns the Eq.(5). We choose the effective area of 
the photon and electron beams intercept to be S e ff ~ S e — 2 • 10~ 4 cm 2 , where 
S e is the transverse area of the LEP2 electron beam, and the total Compton 
cross section a c = 1.86 • 10~ 25 cm 2 at the LEP2 energy. And finally, we should 
take into account the effects of beam mean life time ( we did not consider these 
effects earlier, since every electron bunch provided by the linear accelerator 
could be used only once ). 

In each collision about K N e electrons are scattered. Since every bunch exercise 
the collision with the frequency f pu i se = c/2nR ~ 1.12 • 10 4 Hz, the mean 
lifetime of the beam r b may be estimated as 
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Table 1 

The parameters of the LEP2, SLAC and TESLA electron beams. Here f vu \ se is the 
pulse frequency, H and V are the horizontal and vertical beam radius correspond- 
ingly, iV e is the number of electrons per bunch and is the number of bunches. 





E e (GeV) 


N e (10 10 ) 


n b 


f pulse (Hz) 


H (fim ) 


V ( fim ) 


LEP2 


100 


40 


4 


11.25 • 10 3 


200 


8 


SLAC 


50 


3.5 


1 


120 


2.1 


0.6 


TESLA 


250 


5.15 


800 


10 


0.64 


0.1 



n 



log(l - 8) 1 
log(l - K) f puUe 



400 s, 



where S denotes the maximal fraction of electron loss permitted by beam 
dynamics. In what follows we choose 5 = 0.1, then the number of collisions 
each bunch exercises during one cycle, is equal to 

4.46.10'. 

log(l - K) 

and the mean number of electrons in a bunch N e is given by 



where N e denotes the initial number of electrons in a bunch. 

For the ring accelerators the integrated luminosity per year of operation takes 
the form 



Tint _ -i a7 



n 



r a + n + T f 



frep /3 N e K Tjii 



(8) 



where r a is an acceleration time, t/ is a filling time. Taking into account the 
Eq.(4), one can estimate L % ^ ng as follows 



L -^ 107 -^T7;^^ - 27/rl - 



(9) 
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Table 2 

The real garnma beam parameters at LEP2, SLAC and. TESLA. Here Umax is the 
ratio of the gamma beam maximal energy to the electron energy, # 7 is the opening 
angle of the collimator, (5 denotes the fraction of photons coming through the slit, 
Nry is the total number of converted photons, K is the conversion coefficient and A 
is the laser energy per pulse. 





fj 1 1 I LL Jut 


1 

(10~W) 


P 

1 


l 

(io 4 ) 


K 

(io- 7 ) 


A 

(KT 6 J) 


jj,nt 

(fb- 1 ) 


LEP2 


0.835 


1.26 


0.0265 


0.945 


0.236 


13.4 


0.27 


SLAC 


0.717 


1.93 


0.0202 


1.236 


3.53 


2.88 


0.012 


TESLA 


0.927 


0.759 


0.0422 


0.592 


1.15 


2.21 


0.8 



The electron beam parameters, we have used in the estimations, are contained 
in Tabl. 1. while the 7-beam characteristics are given in Tabl. 2. 

Thus the luminosity at a probable SLAC experiment is about one order of mag- 
nitude smaller compared with LEP2 case due to a lower collision frequency. 
Choosing a more powerful laser one can reach almost the same luminosity 
at SLAC, but the price will be too high: more than one hundred events per 
collision, which reduces significantly the range of the physical phenomena ac- 
cessible for investigations. The luminosity at TESLA even exceeds the LEP2, 
since the TESLA project will operate with multybunch trains. 



3 The Heavy Quark Photoproduction 



The differential cross section of the photon-gluon fusion looks as follows [12] 



da db n _|_ ^ ^ d8a 
di di 9 di 

where X g and £2 denote the gluon and photon helicities. The spin averaged 
and polarized asymmetry distributions in the LO QCD take the form, core- 
spondingly: 
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Fig. 1. The distribution of charm production a C c(pt > Pt) versus p® for differ- 
ent electron energies: E e = 30 GeV - solid line, E e = 50 GeV - dashed line, 
E e = 100 GeV - dotted curve, E e = 250 GeV - dashed-dotted line. The c-quark 
mass is m c = 1.5 GeV/c 2 (a), m c = 1.3 GeV/c 2 (b). 

dcr n naa s (s)e 2 ( 4m 2 s u — m 2 i — m 2 

di s 2 \(i - m?)(u - m?) i - m? u - m 2 

4m 4 s 2 \ , . 

~- 2 2 ■ ( 10 ) 

(t — m 2 ) (u — m?) ) 



dSa naa s (s)e 2 ( ,u — m 2 i — m 2 \2 
di 2 s 2 \\i-m 2 u-m 2 ' 

-s(s - 4 m 2 ) ( 1 + + 1 — )V (11) 

Here s, t, u are the invariant variables of the subprocess, a, a s (s) are the fine 
and strong coupling constants respectively and e q denotes the c-quark charge. 

The produced cc pairs can form J/ip or fragment into D, D*-mesons. In what 
follows we shall consider only the open charm production, because this process 
is more transparent from theoretical point of view. In addition, the open charm 
has an advantage over J / tp production, because its cross-section is at least ten 
times larger for attainable photon energies. 

The accumulated ^-distribution of charmed quarks 



**(P?)= Jd Pt J^^-G(x,Q 2 ) (12) 
versus p® is shown in Fig.l, and Fig. 2 presents the same dependence for 
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0.5 1 1.5 2 2.5 3 3.5 4 
^ p,° (GeV) 



Fig. 2. The distribution of the inclusive D*-meson production <td* x(pt > Pt) versus 
Pt for different electron energies: E e = 30 GeV - solid line, E e = 50 GeV - dashed 
line, E e = 100 GeV - dotted curve, E e = 250 GeV - dashed-dotted line. The c-quark 
mass is m c = 1.5 GeV jc 2 (a), m c = 1.3 GeV jc 2 (b). 

^ 0.6 
^ 0.5 

V s °- 4 

0.3 
0.2 
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-0.2 
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0.1 
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-0.2 
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Fig. 3. The differential asymmetry A^J? (pt) of the c-quark production versus pt for 
the c-quark mass equal to 1.5 GeV/c 2 . The solid (dashed) line corresponds to the 
set A (B) of polarized gluon density [9]. 

inclusive D*-meson production 



In our estimates we have used the parametrizations of both polarized and 
unpolarized gluon densities from [13], c-quark fragmentation functions were 
taken from [14], a s from the global analysis of DIS data [15]. A gamma beam 



30 GeV 


l\ 50 GeV 




5 0.5 1 0.5 1 1.5 2 

p, (GeV) p, (GeV) 


N^^^ GeV 


L 250 GeV 



12 3 2 4 

p, (GeV) p, (GeV) 
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Fig. 4. The differential asymmetry vl^g (Pt) of the c-quark production versus pt for 
the c-quark mass equal to 1.3 GeV/c 2 . The solid (dashed) line corresponds to the 
set A (B) of polarized gluon density [9]. 

spectrum was chosen as monochromatic with energy E 1 = k/(1 + k)E b for 
E e > 50 GeV, while for E e = 30 GeV the smeared-out spectrum, given by 
(2), was used. Numerical integration was performed with the help of adap- 
tive integration code [16]. All distributions have similar form and decrease 
exponentially with p® rise when transverse momentum exceeds 1 GeV j c. One 
concludes that the estimated production rate of D*-mesons with large trans- 
verse momentum pt > 1 GeV jc is sizeable and reaches approximately 30 nb 
(40 nb) at TESLA energy for c-quark mass equal to 1.5 (1.3) GeV jc 2 corre- 
spondingly. 

The differential asymmetry (pt) of c-quark production versus p t 



is shown in Figs. 3-4 for the c-quark mass 1.5 (1.3) GeV/c 2 correspondingly, 
and analogous dependences for D*-meson production are pictured in Figs. 5-6 
for the same mass set. Here we have introduced the following notations 



dp t 
d5a 

dp t 



-I 



ds da n 
s dp t 

ds d5a 
s dp t 



AG(x,Q 2 ). 



As is already known, the differential asymmetry of cc production, derived in 
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Fig. 5. The differential asymmetry AJ^(p t ) of the D* meson production versus pt 
for the c-quark mass equal to 1.5 GeV/c 2 . The solid (dashed) line corresponds to 
the set A (B) of polarized gluon density [9]. 



the LO approximation, has a kinematic zero at p t = m c and then changes 
a sign. For D*-meson production a zero position shifts to the lower value of 
Pt ~ 1 GeV/c due to fragmentation smearing. Integration over total range 
of p t evidently decreases the asymmetry, so it is reasonable to introduce a 
kinematic cut on the D*-meson transverse momentum, say p t > 1 GeV/c. In 
this region AJ^{p t ) does not change sign, moreover the predicted asymmetry 
value heavily depends on the choice of polarized gluon distribution. It is due 
to the fact that the main contribution to production of D*- mesons having 
large transverse momentum comes from the region x > 0.1, where AG(x, Q 2 ) 
is poorely known. At present there are plausible restrictions on the polarized 
gluon density only at x < 0.1, where all model predictions almost coincide. 
Then the precise measurement of polarized gluon density at re > 0.1 will surely 
help to choose a reasonable parametrization of AG(x, Q 2 ). 

Let us briefly discuss the next-to-leading (NLO) corrections to the charm 
production process. As it was shown in [17], these corrections are quite sizeable 
at SLAC and CERN experiments, so one should introduce large K-factor to 
take into account NLO corrections. It is not the case at TESLA, where NLO 
correction to the total asymmetry does not exceed 0.3-0.5. Then one may hope 
that higher order corrections will also be under control at possible TESLA 
experiment. It is not improbable, because all invariants (taking into account 
introduced kinematic cut) are sufficiently large at TESLA energy. The total 
asymmetry is of the order 10 % at \/s Jp ~ 10 GeV and decreases with energy 
rise to a few percent when \fs lp ~ 20 GeV [17]. The introduced kinematic 
cut on D* -meson transverse momentum surely rises this value, because the 
differential asymmetry A]^(p t ) does not change sign beyond p t ~ 1 GeV jc. 
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Fig. 6. The differential asymmetry AJ^(p t ) of the D* meson production versus pt 
for the c-quark mass equal to 1.3 GeV/c 2 . The solid (dashed) line corresponds to 
the set A (B) of polarized gluon density [9]. 

4 Event Reconstruction 



4.1 D*- tagging 

One of the best method for D*- tagging uses the kinematic constraint of the 
decay chain: 



D* — > Dtx — > [Ktx)tx. 
The difference of the invariant masses 



(14) 



AM = m(Kirir) - m{Kn) = m(D*) - m{D) ~ 145 MeV/c 2 



(15) 



is very close to the 7r-meson mass, and can be determined with precision 
about 2.5 MeV/c 2 [6], significantly exceeding accuracy of D*-meson mass 
measurement. It permits to reduce substantially the background to D* + -meson 
production, for example, COMPASS estimates showed that the background is 
less than 10 %. 

The isospin invariance suggests equal D* + and D*° production rates, then 
a(D* + X) ~ 0.5 ■ a(D* X). Since the branching ratio for the decay chain (14) 
is 2.6 %, about ~ 62.4 • 10 4 (83.2 • 10 4 ) D* + -mesons with p t exceeding 1 GeV/c 
will be produced at TESLA energy in a year. Assuming an overall acceptance 
to be equal to 0.2 [6], we expect to reach a reconstruction rate of 
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Table 3 

The target parameters. Here Pt is the target polarization, f t denotes the fraction of 
the polarized target nucleons, T n stands for the target density, and L is the target 
length. 





ft 


Pt 


T n (W 25 cm- 2 ) 


L (cm) 


butanol 


0.24 


0.8 


4 


40 



N D *+ = 124.8 • 10 3 (166.4 • 10 3 ) • year' 1 . (16) 

A rigorous evaluation of acceptance will change these figures, because COM- 
PASS experiment uses kinematic cuts different from ours. Nevertheless, the 
equality (16) gives a good idea about the actual reconstruction rate, attainable 
in the proposed experiment. 

The observed spin asymmetry A^F of _D* + -meson photoproduction with ac- 
count of the beam and target polarizations is given by 



A%> = P b P t f t At. (17) 

Here the asymmetry A]^ is the ratio of helicity dependent and helicity av- 
eraged cross sections for D* + -meson production in 7iV-collisions, calculated 
with account of kinematic cuts. Parameters Pb, Pt are the beam and target 
polarizations, ft denotes the fraction of the polarized target nucleons. The tar- 
get characteristics are given in Tabl. 3, while the gamma beam polarization 
Pb is about unity. 



We impose the p t cut on D*-meson transverse momentum and choose three 
values of p t : 0, 1, 2 GeV/c. The statistical precision of A$ with account of 
target properties is given by 



A<f + ~ = 2.08 (1.81) • 1(T 2 , (18) 

y/l\ D *+ Jt r h f t 

for the cut p t > 1 GeV/c. One can achieve even better statistical precision by 
detecting D* mesons in the total kinematic region of p t (Tabl. 4). On the other 
hand, "strong" cut on the transverse momentum, say p t > 2 GeV / c, will lower 
the statistical accuracy by a factor of three, because the corresponding cross 
section does not exceed 3-4 nb. Our opinion is that the "moderate" restriction 
p t > 1 GeV/ c is best suited for investigation of gluon polarization. Under this 
cut the total cross section is still sizeable, which results in a reasonable value 
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Table 4 

The registration rates of D* + mesons at TESLA. Here a is the cross section of inclu- 
sive D* meson production, N D *+ denotes the total number of observed D* + mesons 
per year with account of efficiency, AAJ^ + stands for the statistical precision. 





a(nb) 


iV D * + (10 3 ) 


A< + (10- 2 ) 


Pt > o 


126 (215) 


524.2 (894.4) 


0.72 (0.55) 


Pt > 1 GeV/c 


30 (40) 


124.8 (166.4) 


1.47 (1.28) 


Pt >2 GeV/c 


4.1 (5) 


17.1 (20.8) 


3.99 (3.61) 



of statistical precision, moreover it permits to decrease background and reduce 
uncertainties in the c-quark fragmentation process. 

4.2 D° tagging 

Experimentally, total number of D°-mesons per charm event is approximately 
equal to the sum of D* + - and _D*°-meson production rates. Produced D° may 
be detected via the simplest two - body decay 

with the branching ratio of 3.8%. An estimate shows that D°- meson recon- 
struction rate is more than corresponding quantity for D* + -meson by a factor 
of three due to a larger value of both the _D°-production cross section and 
decay branching ratio. However, unlike the -D* + -tagging, a background to the 
_D°-production remains significant even after account of kinematic cuts, and 
it exceeds signal, for example at COMPASS, by a factor of four [6]. It means 
that correct evaluation of both the asymmetry and statistical precision of mea- 
surement requires a detailed analysis of the background, which is beyond the 
scope of our paper. 

4.3 Single muon tagging 

The most simple method to select charm production events consists in the 
detecting of muons from semileptonic .D-meson decays. Muons coming from 
light meson decays, as well as Bethe-Heitler pairs, contribute mainly to the 
low values of p t , so one may reduce these backgrounds by imposing cut on 
.D-meson transverse momentum, say p t > 1 GeV / c. Assuming the branching 
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ratio of the decay 



to be 6.8 % for D° and 17.2 % for D + -meson, and ratio of the D-meson 
production rates per charm event to be 



one expect about iV M = 4.51 (6.02) • 10 6 prompt muons at TESLA for m c = 
1.5, (1.3) GeV/c 2 correspondingly. In average, muon acquires a transverse 
momentum equal to about one half of the parent meson mass, that is 1 GeV/ c 
in our case. To obtain a rough estimate, suppose that only one half of muons 
coming from decays of .D-mesons with p t > 1 GeV/c will get a transverse 
momentum larger than 2 GeV jc. In this case a reconstruction rate is given by 



~ 2.26 (3.01) • 10 6 , 

under assumption of 100 % muon detection efficiency The precision of the 
asymmetry measurement with account of the target properties can be roughly 
estimated as 



AArf ~ = 3.5 (3.0) • 10~ 3 . (19) 

The exact evaluation of the statistical error requires a close investigation of the 
background. Nevertheless Eq. (19) gives a good idea about an actual precision 
one may hope to reach at TESLA experiment. We remind, for comparison, 
that SLAC collaboration expected to obtain statistical accuracy of asymmetry 
measurement about 6-10 -3 under severe conditions of large background coming 
from pile up of the events. 



5 Conclusion 



In the present paper we have considered an opportunity to use polarized real 
photon beam for investigation of polarized gluon distributions. Our estimate 
back up the possibility to achieve a higher accuracy in the measurement of 
charm photoproduction asymmetry compared with planned experiments at 
SLAC and CERN. The proposed experiment may give unambiguous informa- 
tion about both the total value and the x- dependence of gluon polarization, 
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which will permit to reduce significantly the number of acceptable models 
decsribing nucleon spin effects. 
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